Spatial Visualization

Motivation

https://research.stlouisfed.org/publications/economic-synopses/2013/08/13/mind-the-regional-output-gap/

trying URL 'http://www2.census.gov/geo/tiger/TIGER2017/STATE/tl_2017_us_state.zip'
Content type 'application/zip' length 9284383 bytes (8.9 MB)
==================================================
downloaded 8.9 MB

number of columns of result is not a multiple of vector length (arg 1)3 parsing failures.
row # A tibble: 3 x 5 col     row   col   expected    actual             file expected   <int> <chr>      <chr>     <chr>            <chr> actual 1  1441  <NA> 57 columns 1 columns 'QGSP_ALL_R.csv' file 2  1442  <NA> 57 columns 2 columns 'QGSP_ALL_R.csv' row 3  1443  <NA> 57 columns 1 columns 'QGSP_ALL_R.csv'

# plot the shapes to ensure we have them
leaflet(states) %>%
  addTiles() %>%
  addPolygons() %>%
  setMapWidgetStyle()

sf_NPR1to1 <- mutate(sf_NPR1to1, FID = as.integer(FID))
sf_USGDP <- dplyr::left_join(sf_NPR1to1, mutate(state_gdp, by = c("FID"="id")))

Error: `by` required, because the data sources have no common variables

# plotting
# ggplot(gdpDF, aes(long,lat, group=id, fill= '2016Q1')) + geom_polygon()
ggplot(state_gdp, aes(map_id=id,fill='2016Q1')) +
  geom_map(map=state_map, color='#ffffff', size=0.1) +
  expand_limits(x=state_map$long,y=state_map$lat) +
  theme_map()

qpal <- colorBin("Blues", domain = sf_USGDP$'2016Q1', bins = 8)
labels <- sprintf(
  "<strong>%s</strong><br/>%g dollars",
  sf_USGDP$state, sf_USGDP$'2016Q1'
) %>% lapply(htmltools::HTML)
leaflet(
  sf_USGDP,
  options= getLeafletOptions(-1.5, -1.5)) %>%
  addPolygons(
    weight=2,
    fillColor= ~qpal('2016Q1'),
    group = 'states',
    fillOpacity = 0.6, opacity = 1, color="#666666",
    highlightOptions = highlightOptions(weight = 4),
    label = labels,
    labelOptions = labelOptions(
      style = list("font-weight" = "normal", padding = "3px 8px"),
      textsize = "15px",
      direction = "auto")) %>%
  addLegend(pal = qpal, values = ~'2016Q1', opacity = 0.7, title = NULL,
    position = "bottomright")%>%
  addLabelOnlyMarkers(
    data=sf_USGDP.centers,
    label = ~as.character(state),
    labelOptions = labelOptions(
    noHide = 'T', textOnly = T,
    offset=c(-4,-10), textsize = '12px')) %>%
  setMapWidgetStyle()

Error in cut.default(x, binsToUse, labels = FALSE, include.lowest = TRUE,  : 
  'x' must be numeric

ggplot(state_map) +
  geom_sf(aes()) +
  scale_fill_viridis() +
  ggtitle("Area of States")

EXAMPLES

# From http://leafletjs.com/examples/choropleth/us-states.js
qpal <- colorQuantile("Blues", domain = state_gdp$'2016Q1', n = 7)
leaflet(states) %>%
  setView(-96, 37.8, 4) %>%
  addPolygons(
    color = ~qpal('2016Q1'),
    opacity = 1,
    dashArray = "3",
    fillOpacity = 0.7,
    smoothFactor = 0.2,
    highlight = highlightOptions(
      weight = 5,
      color = "#666",
      dashArray = "",
      fillOpacity = 0.7,
      bringToFront = TRUE),
    label = labels,
    labelOptions = labelOptions(
      style = list("font-weight" = "normal", padding = "3px 8px"),
      textsize = "15px",
      direction = "auto")) %>%
  addLegend(pal = qpal, values = ~'2016Q1', opacity = 0.7, title = NULL,
    position = "bottomright")

Error in cut.default(x, binsToUse, labels = FALSE, include.lowest = TRUE,  : 
  'x' must be numeric

leaflet(
  sf_NPR1to1,
  options= getLeafletOptions(-1.5, -1.5)) %>%
  addPolygons(
    weight=2,color='#000000', group = 'states',
    fillOpacity = 0.6, opacity = 1, fillColor= ~getFactorPal(state)(state),
    highlightOptions = highlightOptions(weight = 4)) %>%
  addLabelOnlyMarkers(
    data=sf_NPR1to1.centers,
    label = ~as.character(state),
    labelOptions = labelOptions(
    noHide = 'T', textOnly = T,
    offset=c(-4,-10), textsize = '12px')) %>%
  setMapWidgetStyle()